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ABSTRACT 



Many biological electron transfer (ET) reactions are mediated by metal centres in proteins. NADH:ubiquinone 
oxidoreductasc (complex I) contains an intramolecular chain of seven iron-sulphur (FeS) clusters [1], one of the longest 
chains of metal centres in biology and a test case for physical models of intramolecular ET. In biology, intramolecular 
ET is commonly described as a diffusive hopping process, according to the semi-classical theories of Marcus and 
Hopficld 2,3. However, recent studies have raised the possibility that non-trivial quantum mechanical effects play 
a functioning role in certain biomolecular processes [4j [5] . Here, we extend the semi-classical model for biological 
ET to incorporate both semi-classical and coherent quantum phenomena using a quantum master equation based 
on the Holstein Hamiltonian 6J. We test our model on the structurally-defined chain of FeS clusters in complex I. 
By exploring a wide range of realistic parameters we find that, when the energy profile for ET along the chain is 
relatively flat, just a small coherent contribution can provide a robust and significant increase in ET rate (above 
the semi-classical diffusive-hopping rate), even at physiologically-relevant temperatures. Conversely, when the on-site 
energies vary significantly along the chain the coherent contribution is negligible. For complex I, a crucial respiratory 
enzyme that is linked to many neuromuscular and degenerative diseases [7|, our results suggest a new contribution 
towards ensuring that intramolecular ET does not limit the rate of catalysis. For the emerging field of quantum 
biology, our model is intended as a basis for elucidating the general role of coherent ET in biological ET reactions. 



In NADH:ubiquinone oxidoreductase (complex I), the electron transfer (ET) reaction between NADH and quinone 
is mediated by seven iron-sulphur (FeS) clusters that link a flavin mononucleotide to the quinone binding site (see 
Fig. [IJ. Following reduction of the flavin by hydride transfer from NADH, electrons are transferred sequentially along 
the cluster chain, to the quinone. The free energy released by the redox reaction is conserved in transmembrane 
proton translocation, in a separate domain of the enzyme. Thus, ET along the cluster chain must be both fast-enough 
that it does not limit the rate of catalysis, and energetically efficient [8] . Here develop a general treatment of ET in 
metalloproteins, exploiting complex I as a template to elucidate whether quantum mechanical contributions play a 
role biological ET. In our model, each cluster in the chain can exist in two states, oxidised or (one-electron) reduced, 
and each electronic state has a set of vibrational states associated with it. We preserve essential features of the Marcus 
and Hopfield theories [21 E], in which vibrations and vibronic couplings play a central role, by basing our model on 
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FIG. 1: The FeS clusters in complex I connect the NADH/flavin electron donor module (D) to the quinone electron acceptor 
(the electron sink, A). Edge-to-edge distances between the clusters are marked in A, with tunneling amplitudes given in brackets 
in GHz. The structures of the 2Fe and 4Fe cores are inset. The whole assembly is embedded in the protein matrix. Based on 
2FUG.pdb [J. 
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the closely related Holstein Hamiltonian 
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In the first term of Eq. (|TJ) , is the on-site energy of an electron on the i th cluster (2?j does not include interactions 
between the electron and cluster, see below), and c\ and e, are the electron creation and annihilation operators. The 
second term describes the coherent tunneling of the electron between the j th and (j + l) th clusters with amplitude 
tj.j+i (h is Planck's constant). Coherent tunneling gives rise to a continuous wavelike motion of the electron, allowing 
it to exist in a quantum superposition of states across adjacent redox sites. This electron delocalisation is absent 
from diffusive models of ET, such as those of Marcus and Hopfield, in which an electron hops between well-defined 
locations with a certain probability. The third term denotes the vibrational contribution to the energy of the system, 
where Vi is the characteristic frequency of vibration of the i th cluster and the bosonic operators a\ and a* act to raise 
and lower, respectively, the number of excitations in that vibrational mode. The Hamiltonian in Eq. ([I]) describes 
the case of a single vibrational mode per site whereas, in reality, each FeS cluster has a large complement of modes. 
Generalisation of our model to a greater number of modes is possible using the approach outlined in the appendix. 

Associated with the Holstein Hamiltonian is the concept of the polaron, a bosonic quasiparticle formed by the 
electron and the local vibrational modes of the FeS cluster. Polaron formation arises from vibronic coupling, through 
which the reduction of a given cluster leads to a change in its vibrational state; it is described by the fourth term of 
Eq. ([I]) for a vibronic coupling strength of gi at the i th cluster. The polaron concept simplifies both the physical and 
mathematical picture considerably by treating the charge and its accompanying vibrational modes as a single entity. 
The concept is made mathematically explicit by the Lang-Firsov transformation [3] (see appendix), allowing Eq 
to be re-expressed, with renormalised parameters, by Eq. ([2]). 
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Where E* = Ei — Ai is the reduction potential of the i th cluster, Xi = exp ^— (a\ — a ij^j arL d Ai = 

Note that the Hamiltonian in Eq. ([2| is general enough to be applied to a wide range of metalloproteins in which 
ET depends strongly on a set of parameters that can only be extracted from experimental studies. In this study, 
complex I is used as a template to investigate the possibility of quantum mechanical contributions to biological ET 
reactions. Two reduction potential profiles were examined. First, when the seven-cluster chain in complex I contains 
four electrons they occupy alternate positions along the chain |10j . suggesting an alternating reduction potential 
profile, with significant variation in potential between adjacent sites. Alternatively, the same pattern of occupancy 
could result from electrostatic interactions between the sites, suggesting a flatter underlying profile. The relative 
contributions of these two effects remains unclear. For the alternating profile [8], the E* values fluctuate significantly 
(on the order of 0.1 eV) from site-to-site. For such large energetic disorder, the quantum mechanical contribution to 
ET is negligible and the dynamics are well described by purely diffusive modelling. However, using a flat profile in 
which E* is taken to be equal for all clusters other than site 7 (that is known to have the highest reduction potential) , 
the quantum mechanical contributions to ET become significant. Thus, we have chosen this flat profile (with sites 
to 1 to 6 at eV and site 7 at —0.15 eV (36.4 THz)) as the main focus of our study; the flat profile has also been 
adopted in other modelling studies [HI [12] . 

The vibrational modes of protein-bound FeS clusters have been characterised using nuclear resonance vibrational 
spectroscopy, resonance Raman spectroscopy and density functional theory (DFT) [T3] ; we select a common frequency 
vi = v = 334 cm -1 (10 THz), i = 1, 2 ... 7 for all sites in the chain. We estimate the vibronic coupling strength from 
DFT simulations of the so-called 'inner sphere reorganisation energy' A; n by noting the equivalence of this quantity 
to Aj. Ref. [13] provides an inner sphere reorganisation energy for FeS clusters A in = 0.2 eV (48 THz) to yield a 
common vibronic coupling strength for all sites, gi = g = \J v\ m = 22 THz. 

The tunnelling amplitudes are obtained from previous studies of the transfer integrals between FeS clusters within 
complex I using DFT and semi-empirical electronic structure methods [11] . Due to the disorder in spacing between 
clusters, amongst many other factors, these tunneling amplitudes vary between adjacent sites, introducing disorder to 
the Hamiltonian. The tunneling amplitudes all lie within the range tjj+i = 1-95 GHz (see Fig. [I] details), allowing 
the polaron tunnelling term of Eq. ^ to be readily treated as a perturbation. 

The presence of the bosonic operators Xi in Eq. ^ gives rise to two distinct processes that contribute to the 
perturbation expansion [B]. The first of these are so-called 'diagonal transitions' in which the vibrational quantum 
numbers in the chain do not change upon the tunneling of an electron. Diagonal transitions facilitate quantum 
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coherent ET and lead to a modulation of the tunneling amplitude between adjacent sites due to the effect of vibronic 
coupling 

th+i = tjj+iexp ("fa (1 + 2 (M T )») • (3) 

Where (n v (T)) = (exp(his/k B T) - 1) is the thermal average of the number excitations in a vibrational mode with 
frequency v at temperature T and kg is Boltzmann's constant. We thus define an effective Hamiltonian, correct to 
the first order (and assuming a single vibrational mode per cluster), describing the coherent transport of a polaron 
among the FeS clusters 
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The effect of the vibronic coupling has been wholly absorbed into the parameters t* J+1 and E* , resulting in exponential 
suppression of the tunnelling amplitude and a modulation of each on-site energy due to the deformation of the cluster 
by an electron. The second contribution to the perturbation series arises from 'non-diagonal transitions' in which the 
vibrational quantum number may change. The quantum coherence of an electron is not preserved in non-diagonal 
transitions, allowing them to be described stochastically by a set of hopping rates that lead to diffusive motion. These 
rates are calculated within our model using the Marcus- Jortner equation |14j 
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and by enforcing detailed balance for the reverse rates. Here, v =10 THz is a typical frequency for the optical 
vibrational modes of the environment and A ou t is the outer sphere reorganisation energy describing the interaction of 
the polaron with its surrounding protein scaffold and solvent. There is no consensus value for A ou t for FeS clusters 
in the existing literature. We assume that the outer and inner sphere reorganisation energies are of comparable 
importance by choosing A ou t = Aj n = 0.2 eV (48 THz). 

The incoherent hopping rates Eq. ^ enter into the description of ET within the chain as Liouvillian terms within 
a Born-Markov quantum master equation alongside the Hamiltonian dynamics (Eq. (EJ)) (see appendix for details). 
Further terms are included in the master equation to address the fact that the closed system we have considered thus 
far is an inadequate description of 'in situ' ET. In complex I, the FeS clusters are embedded within the protein scaffold 
and surrounding solution that constitutes their environment. Interaction of a quantum system with its environment 
introduces non-reversible processes to its dynamics; these manifest themselves as local dephasing - a process which 
reduces quantum superposition states to classical statistical mixtures by the process of 'decoherence' at a given rate. 
In addition, we include the irreversible dissipation of the polaron from the final FeS cluster of the chain to the quinone 
at a rate R = 1 GHz that is of the order of other ET processes in the chain. The origin of local dephasing can be 
associated with the outer sphere reorganisation energy A ou t . In lieu of a detailed knowledge of the density of modes 
for the environment, we assume a dephasing rate of the form 

7 = 7 o(l + K a (T)))« 7o (l + ^, (6) 

where 2/ a = 0.1 THz is taken as a typical frequency for the vibrational modes of the environment that are responsible for 
dephasing |15j and 70 is the projected magnitude of quantum dephasing at zero temperature. Considering Eqs. ([3]), ^ 
and @, we see that with increasing temperature the amplitude of coherent tunneling is reduced as the rate of 
incoherent tunneling and local dephasing increases. 

The time evolution of the system is obtained by standard numerical integration of the Born-Markov master equation 
with the boundary condition that the electron is initially localized at the first site of the chain. The motion of the 
electron among the seven site chain is governed by a number of factors, notably; the disorder of the chain - which here 
enters via the tunneling amplitudes - and the temperature of the system. At low temperature, the coherent (diagonal) 
contribution to ET is significant compared to the rate of incoherent (non-diagonal) processes. As the temperature is 
increased, the importance of coherent transport is reduced due to the suppression of the coherent tunneling amplitudes 
and the onset of dephasing and incoherent tunneling, leading to diffusion-like transport. For high enough values of 
the zero temperature dephasing (70 ~ 1 GHz) at T — 300 K, quantum coherence is no longer maintained on a 
timescale relevant to ET and purely diffusive transport is recovered with rates given by the Marcus- Jortner equation 
(see appendix). Fig. [2] shows the sensitivity of the transport to the temperature and amount of zero temperature 
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FIG. 2: The time taken to achieve 50 % probability of quinone reduction (50% population of the sink) at 310 K, after the 
electron was initialised in site 1, as a function of zero temperature dephasing 70. The shortest time for quinone reduction, 
35 /xs, is observed at 70 = 0.004 GHz. This is almost 10 times shorter than the analogous semi-classical treatment using the 
Marcus- Jortner equation (304 /is). For comparison, diffusive modeling using the Marcus equation [2] predicts a time of 136 /is 
and the empirical treatment of Dutton et al. [12] estimates the time required to be 7 (j,s (for a reorganisation energy of 0.4 eV). 
Inset: The probability of quinone reduction is shown at t = 40 /j,s as a function of 70 for various temperatures. 

dephasing 70. 

At 10 K there is little probability for the electron to reduce the quinone with any level of dephasing. This is a 
feature of the reduction potential profile of the chain: While sites 1-6 have similar on-site energies, the energy of site 
7 is 0.15 eV (36.4 THz) lower. This difference in energy is much greater than the tunneling amplitude between sites 6 
and 7. Consequently, coherent transport is strongly suppressed between these sites and occupation of site 7 is achieved 
by non-diagonal terms that become more pronounced at higher temperatures. Using this profile, the final step in the 
ET chain of complex I is thus a thermally activated incoherent hop, though transport elsewhere in the chain may be 
quantum coherent. The alternating profile discussed above yields purely incoherent dynamics for similar reasons. 

The rate of ET from site 1 of the chain to the quinone molecule increases between 10 - 150 K across the full range 
of 70 that was investigated. Above 200 K the transport becomes increasingly sensitive to the value of 70- At 310 K 
the fastest electron transfer to the quinone is achieved for 70 ~ 0.004 GHz (see Figs. [2]-|4| but for other values of 
70 the transport is far less efficient. These features are akin to the 'dephasing assisted transport' 16 18J that has 
been much discussed in the context of quantum coherent exciton transport in photosynthetic complexes. Dephasing 
assisted transport predicts that an intermediate region between purely coherent and incoherent dynamics offers the 
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FIG. 3: The populations of the seven sites of the FeS cluster chain in complex I at T = 310 K and 70 = 0.001 GHz. The timescale 
for ET to the sink (reduction of the quinone) is roughly 10-100 fj,s, over which the ET appears diffusive. The populations of 
sites 1 & 2, and 3 & 4, are overlaid, while sites 5-7 are barely ever populated. Much shorter timescales (100 ns) show the 
presence of quantum coherent transport between the first two sites of the chain (inset). For comparison, the population of the 
sink is shown for the case of purely diffusive dynamics calculated using the Marcus- Jortner equation (dashed line). 

fastest rate of transport. Notably, the behaviour may not be solely attributed to dephasing, due to the presence 
of incoherent tunneling, but the notion of ET enhancement by the interplay of quantum coherent and incoherent 
processes remains. 

The presence of dephasing assisted transport, and hence partly quantum coherent ET within complex I, at phys- 
iological temperature and in certain regimes of the parameter 70, is further explored in Fig. |3j which shows the 
population of each of the 7 sites of complex I as a function of time for 70 ~ 0.001 GHz at 310 K (similar plots for 
different values of 70 are given in the appendix). On the timescale of ET from site 1 to the quinone molecule (10 - 100 
fis) the transfer appears qualitatively diffusive. The electron quickly equilibrates among the first four sites - which are 
to a certain extent decoupled from the rest of the chain due to the large distance between sites 4 and 5 - by thermally 
activated hopping. Upon making a transition from site 4 to site 5, the electron is more likely to proceed toward the 
sink than to undergo a transition back to site 4. The disorder in the tunneling amplitudes therefore introduces an 
effective directionality to the electron motion. Further directionality arises from the thermally activated hop between 
sites 6 and 7, and the low population in sites 5-7 may be described by the electron quickly being transferred to the 
sink upon transition from site 4 to 5. On much shorter timescales (100 ns), coherent oscillations of the electron persists 
between sites 1 and 2 of the chain (see inset to Fig. [3]) and to a lesser extent between sites 3 and 4. These oscillations 
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FIG. 4: The population of the first four sites of the cluster chain and the quinone as a function of temperature for 70 = (solid 
lines) and 'optimal' dephasing (70 = 0.004 GHz, dotted lines) after 40 /is. 



are the hallmark of quantum coherence and are contrary to the typical relaxation-like behaviour exhibited in classical 
transport. Within our model, coherent oscillations can persist for upwards of 100 ns for 70 < 0.004 GHz. Remarkably, 
despite quantum coherence only persisting on a timescale two orders of magnitude slower than that of the overall ET, 
quantum effects have a significant effect on the transport rate, yielding roughly ten times the probability of quinone 
reduction after 40 /is compared to the analogous purely diffusive Marcus- Jortner model. To further emphasise this 
feature, Fig. [4] shows the population of the quinone molecule and sites 1-4 (as sites 5-7 are only ever sparsely 
populated) when the chain is subject to different rates of dephasing. For no dephasing (70 = 0), the probability of 
quinone reduction does not exceed 0.25 within 40 /is for T < 310 K. However, for 70 = 0.004 GHz the rate of transfer 
is increased at all temperatures, with a maximum probability of roughly 0.7 for quinone reduction after 40 /is. This 
is evidence that the addition of local dephasing processes can enhance ET that is already partly incoherent due to 
non-diagonal processes. 

We have presented a full microscopic quantum model of ET in metalloproteins using a novel approach based on 
the Holstein Hamiltonian [5] and a quantum master equation. By accounting for the unique properties of quantum 
dynamics up to the first order of perturbation theory, we show that, in principle, even small quantum mechanical 
contributions can enhance ET rates in metalloproteins by almost 10 times, compared to analogous semi-classical 
models under certain conditions. Our model also extends the study of ET to temperatures that are too low for 
thermally activated stochastic hopping models to apply. For a flat reduction potential profile in complex I, we 
predict an increased ET rate at moderate dephasing that is attributable to quantum coherence and its interplay with 
decoherence processes. For reduction potential profiles that vary markedly from site-to-site, our model converges with 
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well-established semi-classical treatments of ET [T3]. We note that reduction potential measurements take place on 
far longer timescales than ET. Consequently, the 'transient' reduction profile relevant to ET may not include energetic 
contributions from slow relaxation processes. This could reconcile the requirement of quantum coherence for a flat 
reduction potential profile with the disorder in reduction potentials generally observed in biological molecules. The 
quantum coherent enhancement of the ET rate in complex I is reminiscent, though not directly analogous to, recent 
studies of coherent exciton dynamics within photosynthctic complexes |16H19j . Our model is general enough to be 
applied to other metalloproteins and can be adapted to investigate other non-biological ET processes also. 
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Appendix A: Lang-Firsov transformation 

We show that a more general Hamiltonian may be used as the starting point of our study than that of Eq. 
We also note that it may serve as the starting point to investigate the ET dynamics of metalloproteins other than 
NADH: ubiquinone oxidoreductase (complex I). The extended Holstein Hamiltonian [5] can be partitioned into elec- 
tronic, vibrational and interaction components (H c \, H V ib and H m t respectively) as 

H = H cl + H vih + H int . (Al) 

The electronic part of the Hamiltonian describes an N site tight-binding lattice with disorder in both the on-site 
energies and the nearest-neighbour tunneling amplitude ij.j+i 

N N-l 
* 3 

The vibrational component, flvib, describes a set of harmonic oscillators associated with each of the N lattice sites. The 
mode at site i has a characteristic frequency Vi : ki and corresponding bosonic creation and annihilation operators, 
a\ k . and a^fc. 

N 
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Finally, the interaction term H m t accounts for the vibronic coupling 
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Where gi.ki is the coupling strength between the kf 1 phonon mode and an electron at site i. The vibronic coupling 
thus describes the displacement of the phonon modes at a given site conditioned on the presence of an electron at 
that location. 

Combining Eqs. (All - (A4) we obtain the full Holstein Hamiltonian 
H = H c \ + H vi b + Hi D t ■ 
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The Holstein Hamiltonian is often treated using a polaron approach, facilitated by the Lang-Firsov [9] (or 'polaron') 
transformation using the unitary operator 

V = e s , 

^EE^ti-^)' ( A6 ) 



The transformed fermionic operators are calculated with use of the Hadamard lemma 1 

c } = Vc 3 V\ 
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Where Eq. (A8| follows from hermitian conjugation of Eq. (A7). 
Similarly, the transformed bosonic operators are 
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With Eq. (AlOl following from hermitian conjugation of Eq. (A9|. Combining Eqs. (A7| - (AlOl, the Lang-Firsov 
transformed Holstcin Hamiltonian is found to be 

H = VHV\ 
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Appendix B: Diagonal processes 

We provide a derivation of the renormalised coherent tunneling amplitude arising from 'diagonal processes' - in 
which the number of phonons remains unchanged within the chain upon a tunneling of an electron. The following 
derivation follows that of Ref. [20] closely. Beginning with Eq. ^ 

N N-l N 
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Where, now 

ha 2 

A (Bl) 
Xi = a® (^-^ (4 - a^} . (B2) 



For A, B G C mxm , the Hadamard lemma states: e A Be~ A = B + [A, B] + i [A, [A, B]] + A [A, [A, [A, B]]] + ... 
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Considering the case when t/j+i <C -Ej, z^, the polaron tunneling term of Eq. ^ may be treated as a small perturba- 
tion. To first order of the perturbation expansion we must evaluate terms of the form 
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Where \{n}) (\{m})) is a multi-mode Fock state describing the phonon occupation numbers in the initial (final) state 
by the set of numbers {n} = rii, n^.-.n^ ({ m } = m i, ™2--- m Jv)j aud |z) denotes the state of an electron located at site 
i. The diagonal contribution t* i+1 to the object in Eq. (B3) amounts to the evaluation of the following quantity 
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Factoring the fcrmionic and bosonic parts of Eq. ( B4 ) gives 
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Where the second line follows from the orthonormality condition 



Haussdorff-Zassenhaus formula 2 to Eq. (B5) gives 
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Where, the \rii) are single mode Fock states describing the vibrational state of the i th site. To proceed, the thermal 
average of the quantity in Eq. ( B6 ) is taken, 
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is the partition function for the i vibrational mode 



1 oo 



n— v 



(n|exp (£a^) exp (— £a) |n), 



It can be shown (see Ref. [20] for further details) that 
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where (n(T)) = (exp (pf ) — l) 1 is the Planck average of the number of phonons. Substituting this result into 
Eq. (B7), finally gives the thermally averaged rate of diagonal transitions 
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2 For A, S G C mxm , r, t 6 C, [A, B] = r, the BCHZ formula reduces to: exp (t (A + B)) = exp (tA) exp (tB) exp (-■ 
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We see that the vibronic coupling manilests itself as an exponential suppression of the polaronic tunneling amplitude 
tjj+i relative to its 'bare' electronic counterpart £j - ,j+i. This is attributable to the effective mass acquired by the 
electron due to its coupling to phonons in the polaron picture. Polaron tunneling is also exponentially sensitive to the 
number of phonons excited in the relevant clusters; a quantity that increases as a function of temperature, leading to 
further suppression of the tunneling amplitude t* 



Appendix C: Born-Markov master equation 



Interaction of a quantum system with an environment introduces non-unitary processes to its dynamics. These 
manifest themselves in the form of local dephasing processes; that lead to the 'decay' of superposition states to 
statistical mixtures of those states, and dissipation; the irreversible loss of a particle from the system. The description 
of these 'incoherent' processes is aided by adopting the density matrix formalism 3 . This is convenient when describing 
dissipative processes, as the diagonal elements of the density matrix can be identified as the 'population' of the sites, 
while the off diagonal elements describe the coherence between different basis states. 

The origin of local dephasing within our model can be associated with the outer sphere reorganisation energy 
A out describing the interaction of the polaron with its surrounding protein scaffold and solvent. We argue that 
environmental vibrational modes interact only weakly with the polaron due to the screening of the electron charge by 
the distorted cluster ions, allowing us to invoke the Born approximation. The additional assumption of a Markovian 
environment allows the dephasing process to be written in the form of a Lindblad term 

N 

£dc P h = Z)7 (- { c Kp(*)} + 2cjcip(i)cjci) . ( C1 ) 

i=l 

Where 7 = 70 (1 + ksT/hv) is the rate of dephasing at all sites. In lieu of a detailed knowledge of the density of 
modes for the environment, we treat the projected magnitude of quantum dephasing at zero temperature, 70 as a free 
parameter and scan the range 70 = — 1 GHz to investigate the robustness of the polaron transport under dephasing. 

The ultimate goal of ET in complex I is the reduction of a quinone molecule located at the end of the Fe-S chain. 
Once the quinone is reduced it dissociates from the enzyme. In this way the quinone acts as a sink and introduces a 
further non-reversible process to the dynamics of the system. Dissipation of the polaron from the N th site at a rate 
R is mediated by another Lindblad term 
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(p(t)) = R (- {4c 7 , p(t)} + 2c 7P (t)4) . (C2) 



Dissipation from all other sites can be accounted for in a similar way, however the timescale for this process is taken 
to be much longer than the transport time of the electron, and as such is irrelevant to our model. 

Incoherent tunneling enters our model via the Marcus- Jortner rate equation (see Eq. ([5])). This process accounts for 
transitions of the electron between two sites in which its phase coherence is not preserved, due to inelastic scattering 
with phonons. For sites j and j + 1, the Marcus- Jortner equation predicts a forwards rate Oj-^+x and, by enforcing 
detailed balance, a backwards rate Oj+x-y; of electron transfer. The Liouvillian term describing incoherent tunneling 
is thus 

N-l 

Aunn (p(*)) = a J->i+l ( _ { C ] C J+lcJ +1 Cj-p(t), /?(*)} + 2c ] + l C ^( i )4 C i + 1 ) 

3=1 

+ aj+i^j (- { c ]+iC]c]c j+1 p(t), p{t)^ + 2c)c j+l p{t)c) +l c^j . (C3) 
The full evolution of the density matrix is then described by the Born-Markov master equation 

^ = - \ [H eS , P {t)\ + C dcph (p) + £ sink (p) + C tunn Go) . (C4) 
Where H c ff is the Hamiltonian given in Eq. @. 



3 The reduced state of the system refers to the density matrix obtained after a partial trace is performed over the environmental degrees 
of freedom for the joint denisty matrix describing the state of the system and environment together (see for example I21p . In the case 
of complex I the state of the system may be written as |V>(i)) = ^j— bj(t)\j), where the bj(t) are time dependent amplitudes. The 
reduced density matrix describing the system is then p(t) = \ip(t))(ip(t)\. 
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FIG. 5: The populations of the seven sites of complex I, and the probability of quinone reduction as a function of time with 
zero dephasing at T — 310 K. Coherent oscillations between sites 1 and 2 persist for several microseconds before the onset of 
qualitatively diffusion-like behaviour. The probability of quinone reduction does not exceed 0.25 after 40/is. 

Appendix D: Additional results 

Fig. [3]and Figs.[5]-[6]show the effect of dephasing upon ET within complex I. For a low rate of dephasing (70 < 0.001 
GHz) at physiological temperature (see Fig. [5]), the presence of coherent oscillations between the populations of 
adjacent states is observable on microsecond timescales. These oscillations arc the hallmark of quantum transport. 
The amplitude of the oscillations is damped primarily by the onset of incoherent tunneling and ET qualitatively 
approaches diffusive behavior on the longer timescale of quinone reduction. For high rates of dephasing (70 > 1 GHz 
) and T = 310 K, quantum coherence does not persist on a timescale relevant to ET and the dynamics are purely 
diffusive. Fig. [6] shows that in the limit of large dephasing, the dynamics quantitatively approach the behaviour 
expected using Marcus- Jortner rates within a classical master equation. 

ET at intermediate rates of dephasing (70 ~ 0.001 GHz) and T = 310 K exhibit coherent oscillations that persist 
for approximately 100 ns (see Fig. [3|, approximately 100 — 1000 times shorter than the typical timescale of quinone 
reduction. There are two noteworthy aspects of the ET in this regime: i) The rate of quinone reduction is significantly 
enhanced compared to the cases of low (more 'quantum') and high (more 'classical') dephasing. ii) Quantum coherence 
plays a significant role in ET despite persisting for only a fraction of the typical timescale of quinone reduction 
(~ 10 — 100/is). The enhancement of transport rate in the presence of intermediate dephasing arises from the 
interplay of quantum coherent and classical transport. This behavior is reminiscent, though not directly analogous 
to, the notion of 'dephasing assisted transport'. In our model, the interplay of classical and quantum transport is 
richer due to the presence of incoherent hopping alongside pure dephasing. 
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FIG. 6: Left: The populations of the seven sites of complex I, and the probability of quinone reduction as a function of time 
under the influence of high dephasing (70 > 1 GHz) at T — 310 K. Quantum coherence is not observed on a timescale relevant 
to ET and the typical Marcus- Jortner diffusive transport (right) is recovered exactly. The probability of quinone reduction 
does not exceed 0.1 after 40/is 
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